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ABSTRACT 

We present analyses of two Ginga observations and two observations from the ROSAT database o 
NGC 7469, focusing on the spectral variability observed on timescales of days and longer. During the . 
day 1988 Ginga observation, the hardness ratio (8-21 keV/3.4-5.7 keV) increased significantly as th* 
total flux decreased by 30%. As the spectrum is well fit by the reflection model and since the spectra 
variability dominates the higher energy band, this could be explained by either a variation in the power 
law index or in the effective covering fraction of the reflecting material. This ambiguity is inherent ii 
reflection modeling of Ginga spectra from moderate flux Seyfert 1 galaxies. Assuming that the power4av 
index did not change, we find that the reflected flux is consistent with being constant, suggesting tha 
much of the reflecting material may be located more than 3 light-days from the continuum source will 
the molecular torus being a plausible site. This scenario is also supported by the report of a narrov 
rather than broad iron Ka line in the ASCA data by Guainazzi et al. NGC 7469 was faint during th^ 
1989 Ginga observation, but variability was observed with doubling timescale of 5 hr, and the spectrun 
was harder. A reflection component could not be constrained, and the change in the spectrum could b 
explained by an increase in neutral absorption. 

The brighter of two ROSAT spectra was significantly softer, and in both spectra there was evidence c 
spectral complexity, as has been previously reported by Turner, George, & Mushotzky and Brandt et al 
The spectrum could be fit by a variety of two-component models, including a warm absorber model, ai 
ionized disk model, and a thermal model with single-component blackbody spectrum, but joint fitting c 
the 1988 average Ginga spectrum and the nonsimultaneous ROSAT spectra favored thermal models, a 
other models required an anomalously high reflection ratio. This model is supported by the observatio: 
of a soft excess component and the lack of ionized absorption edges in the ASCA spectrum by Guain 
azzi et al. The long-term spectral variability could be explained by relative variability between ih 
power-law and soft excess component normalizations, perhaps implying that hard X-ray reprocessing i: 
thermal material does not dominate on long timescales. 

Subject headings: galaxies: individual (NGC 7469) — galaxies: Seyfert — X-rays: galaxies 


1. INTRODUCTION 

NGC 7469 is a moderate-luminosity (La-io ~ 3 x 10^^ 
ergs s“^) Seyfert 1 galaxy, first detected in X-rays using 
Uhuru (Forman et al. 1978). Subsequent observations dis- 
covered rapid, large-amplitude flux variability including an 
observation by Ariel 5 of a 2 day flare in which the flux 
increased by a factor of 6 (Marshall, Warwick, & Pounds 
1981). Rapid variability by a factor of 2 in 5 hr was observed 
in the 1-7 keV X-rays with no accompanying variability of 
the 0.05-2 keV X-rays during one EXOSAT observation, 
while variability in both bands occurred among four 
observations in 10 days (Barr 1986). 

Evidence for a soft excess component has been found 
often. When the EXOSAT spectra, were fitted with an 
absorbed power-law model, the absorption column was 
found to be below the Galactic value, implying the presence 
of a soft excess component (Barr 1986). A soft excess com- 
ponent was also found in the Einstein IPC -h MPC spectra 

* Current address: RIKEN, Institute of Physical and Chemical 
Research, Hirosawa 2-1, Wakoshi, Saitama 351-01, Japan. 


(Urry et al. 1989) and SSS -h MPC spectra (Tur 
1991). ROSAT spectra showed evidence for sof 
complexity (Turner, George, & Mushotzky 1993. 
et al. 1993). A soft excess component was also fo\ 
1993 ASCA spectrum (Guainazzi et al. 1994). Cc 
all the EXOSAT data, Grandi et al. (1992) founc 
spectrum became flatter as the flux increased. 
Courvoisier (1992) also analyzed the EXOSAT sp 
found that they were consistent with a constant ' 
component. Evidence for hard X-ray reprocessing 
found including a hard reflection component an 
line in the 1988 averaged Ginga spectrum (Piro, 

& Malsuoka 1990), w'hich was found to be narr 
ASCA spectra (Guainazzi et al. 1994). Evidence ft 
line was found in the EXOSAT Spectral Survey : 
the 97.5% confidence level (Leighly, Pounds, c 
1989). Strong evidence for a reflection compo 
also found in the HEAO 1 A-2 data (Weaver, / 
Mushotzky 1995). 

In this paper we investigate the spectral variabi 
ROSAT and Ginga data. In § 2 the data, reduc 
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liming analyses arc described including hardness ratios, 
which clearly demonstrate the presence of spectral variabil- 
ity. Section 3 gives the spectral fitting results which are used 
to determine which component likely varies producing the 
observed spectral variability. In § 4 these results arc dis- 
cussed in terms of current physical models of active galactic 
nuclei (AGN). 

2. DATA 

2.1. The Ginga Data 

The Ginga mission and Large Area Counter (LAC) are 
described in Turner et al. (1989), Conservative selection cri- 
teria were used (COR > 9, counting rate above the upper 
discriminator less than 5.95 counts per detector, and earth 
elevation angle greater than 6*^). The shape of background 
spectrum, particularly in low energies, varies as a function 
of the ascending node of the satellite due to radioactive 
decay (e.g., Hayashida et al. 1989), and this can result in 
spurious spectral variability. We employed additional con- 
servative selection criteria in order to minimize possible 
systematic errors. The data were selected further so that the 
same range (within SOC^-ISO®) of ascending node was rep- 
resented in the source and background data. Furthermore, 
spectral variability was investigated for timescales no 
shorter than 1 day, over which the effect of the main com- 
ponent of the variable background shape is expected to 
average out. Remaining background and source data were 
sorted and subtracted according to the count rate above the 
upper discriminator. Only data from the top layer of the 
LAC detector were used for spectral fitting because of 
general overall better statistics. Details of Ginga data selec- 
tion and this background subtraction method, commonly 
used at Nagoya University and RIKEN in Japan, can be 
found in the Appendix of Leighly et al. (1994). Further dis- 
cussion of the reliability of this method is given in Mat- 
suoka et al. (1990), and a discussion of the effect of this 
background subtraction method on spectral variability 
investigations is given in Fiore et al. (1992). 



Fig. 

ratio light curves, 
note the decrease in flux during the last part 
conesponding to transmission T > 0.6. 


NGC 7469 was first observed for 3 days in 1988 from July 
10 through July 13, and source-free sky was observed for 2 
days during 1988 July 8 and 9. The total flux light curve 
(1.1-20.9 keV) is shown in the top panel of Figure la, where 
the zero reference is 1988 July 10 1:11:59.509. During the 
first day, the average flux (top plus midlayer) was about 24 
counts s"S but it decreased after 2.2 days to an average of 
about 18 counts s"^ 

Spectral variability occurred during this observation. The 
hardness ratio computed from the top plus midlayer data 
(8.0-20.9/3.4-5.7 keV), shown in the second panel of Figure 
la, demonstrates that the hardness is variable on the time- 
scale of days with the spectrum hardening as the flux 
decreased, opposite the EXOSAT result (Grandi et al. 
1992). The softness ratio (1.1-3.4/3.4-5.7 keV) varies also, 
but with much less significance. Comparing the average 
ratio for the first day with the average ratio for the second 
and third days shows a reduced for the constant model of 
2.3 for the softness but 10 for the hardness ratios. 

A second observation was made on 1989 November 
23-24, with nearby blank sky being observed on 1989 
November 22. The light curve from the 1989 observation is 
shown in Figure \b, where the zero reference point is 1989 
November 23 19:11:6. Satellite attitude problems occurred 
during this observation, with the transmission decreasing 
from ^0.78 from the beginning of each day, so that before 
aspect correction the average flux was about 5 counts s ^ 
For spectral analysis, only data with collimator transmis 
sion >0.6 were used as marked in Figure lb. During the 
second day, the flux decreased by about 40% in 20,000 s 
This decrease is consistent with the factor of 2 decrease ir 
5-6 hr observed using EXOSAT (Barr 1986). The hardnes 
ratio in the lower panel shows that no spectral variabilit} 
occurred during this observation, but the spectrum is sig 
nificantly harder than was observed in 1988. 

2.2. The ROSAT Data 

Two on-axis ROSAT observations of NGC 7469 wer* 
available in the archive. The source spectrum extractioi 
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Fig. lb 


^ chnu/e 1 1-209 keV total flux and the bottom panel shows the 8.00-20.9 kcV/3.3-5.7 keV hardn 

' '"‘'“reJeSy’ (^Gin^a 1988 observation; note the general trend of increasing hardness wth d^reasing flux, (fc) Ginga 1989 observati. 

^ r t ; tf of the second day of observation. The marked intervals show the data accumulated for spectral anal) 
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Fig. Z—ROSAT light curves, where the top panels show the 01-2 keV total da^'^tuh^th^^ghM Am 1991 

light curves, respectively. Significant fiux variability but no significant ^ 

observation showing a larger softness ratio, (a) ROSj4T 1991 observation, (h) 1992 


radius used was 312, and background was taken from an 
annulus. The “gain2” detector response matrix was used 
for spectral fitting as recommended for observations 
occurring after 1991 October. We are aware of the remain- 
ing calibration uncertainties of the ROSAT PSPC which 
are discussed in detail by Brinkmann et al. (1994) and Fiore 
et al. (1994). The possible effect of this residual uncertainty 
on the observed spectral variability is discussed in § 3.2. 

NGC 7469 was observed on 1991 November 28 for 4251 
s, and the mean flux was 3.0 counts s" *. The light curve in 
Figure 2a shows that between the two observation intervals 
the source flux increased by ~25%. Spectral fitting results 
from this observation were reported previously by Turner et 
al. (1993a), and differences between those results and these 
are probably attributable to the use of the updated detector 
response matrix. NGC 7469 was also observed using 
ROSAT from 1992 May 28 to 1992 May 31 for 18,575_s. 
The average flux from this observation was 1.72 counts s . 
The light curve shown in Figure 2b demonstrates that a 
40% decrease in flux occurred during the observation. Spec- 
tral fitting results from this observation were reported by 
Brandt et al. (1993). The softness ratios (0. 1-0.4 keV/0.5-2.0 
keV) show that there is no significant spectral variability 
during either observation. However, the pulse-height 



Energy (keV) 

YiG. 3. PHA ratios of ROSAT spectra showing that the spectrum 

from the 1991 observation is substantially softer than that from the 1992 
observation. 


analyzer (PHA) ratio given in Figure 3 shows that the spec- 
trum hardened, particularly below 0.5 keV, as the flux 
decreased between the two observations spaced 6 months 
apart. 

3. SPECTRAL FITTING RESULTS 
3.1. The Ginga Data 

The 1988 data were accumulated into four spectra includ- 
ing the total spectrum as well as spectra from each day 
separately, while the 1989 data were accumulated into one 
spectrum. These spectra were fitted from 1.11 to 20.9 keV 
with an absorbed power law plus Gaussian emission-line 
model, with the line energy fixed to 6.4 keV and width (a) 
fixed to 0.05 keV. The parameters from these fits are listed 
in Table 1, with the line flux in column (4) and the line 
equivalent width in column (5). All fits were acceptable, 
except for G88.1, where the poor fit can be partially attrib- 
uted to the lowest energy channels. Ignoring an additional 
lowest two channels results in an acceptable fit (xC = 1-25), 
with no significant difference in continuum fit parameters. 
The absorbing column densities are generally consistent 
with the Galactic value; note that columns below 10 
cm”^ cannot be determined using Ginga. The apparent 
photon index is variable and clearly correlated with the flux, 
as suggested by the hardness ratios shown in Figures la and 
lb. The line flux does not vary significantly, although there 
is the suggestion of a decrease in line flux in the lower flux 
spectrum from 1989. The equivalent width appears to 
increase, but the change is not significant. 

3 . 1 . 1 . Reflection Model Fitting 

The reflection model has been shown to describe Ginga 
spectra from Seyfert 1 galaxies well (e.g., Nandra & Pounds 
1994). Two scenarios represent the extremes of variability 
behavior in the reflection model, and these can be used to 
place constraints on observed spectral variability (e.g., 
George et al. 1993). If the reflecting material is either far 
from the central engine or extends far from the central 
engine compared with the variability timescale, the 
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TABLE 1 

Parameterizing Model Fit Results 


N„ Photon Flu* 

Spectrum (10^‘ cm'^) Index PI* P2'’ (10" ergs cm'* s’*)® 

(1) (2) (3) (4) (5) (6) (7) 


G88.tot <0.3 52lH 130!|g 4.1 30/29 

G88.1 <0.8 1.85!g;g^ S.llU HOl^g 4.8 47/29 

G88 <0.4 1.77!gg5 S.2111 140!jg 3.5 21/29 

G88.3 <0.7 1.77!g;g5 5.2! 130!Jg 3.8 22/29 

G89.2 0.4 <3.3 1.58!g;i? 3.9!^J 180!lJg 2.0 25/29 

R91 0.67!g;gt 2.4 ± 0.1 0.85!g;?J 0.27 ±0.13 3.4 22/25 

R92 0.57 ±0.2 2.2 ± 0.1 0.92!g;gg 0.28 ± 0.09 2.0 11/25 


■ Additional model parameter. For the Ginga spectra, this is the narrow (0.05 keV) line flux (x 10" * photons 
cm“* s" *) at fixed energy 6.4 keV, and for the ROS.4T spectra it is the energy of the absorption edge (kcV). 

** Additional model parameter. For the Ginga spectra, this is the equivalent width of the 6.4 keV iron line 
measured in cV, and for the ROSy4 7 spectra it is the optical depth of the edge. 

* For Ginga, this is the 2-10 keV flux, while for ROSAT it is the 0.1-2 keV flux. 


response of the reflecting material would be lagged and 
smeared, and this would be observed as a constant normal- 
ization of the reflected component. The following model ex- 
presses the reflection component normalization explicitly: 

F{E) = -h A{E, T) + 7^(£)] . (1) 

A{E, T) is the reflection spectrum for a face-on orientation, 
and Ii(E) is the Gaussian expression for the iron Ka fluores- 
cence line, with energy and <t fixed as before. In contrast, if 
the reflection region is small compared with the source 
variability timescale, the ratio between the normalization of 
the reflection component and the power-law component 
referred to as the reflection ratio r should be constant. The 
following model expresses the reflection ratio explicitly: 

F{E) = -h M(£, T)] + I^[E)} . (2) 

A face-on orientation with the reflecting material sub- 
tending a solid angle il — 2n to the primary hard X-ray 
emission is implied when r == 1. 


The fit results to both models are given in Table 2. Piro et 
al. (1990) previously analyzed the average spectrum from 
the 1988 observation; our results are consistent with theirs, 
and the smaller fit xl horn this analysis probably reflects a 
more conservative data selection. The low flux of the 1989 
observation resulted in poor statistics, so that the param- 
eters of the reflection model could not be constrained. 
Overall, the values obtained for the reflection ratio are gen- 
erally larger than 1, though most are consistent with 1 at the 
90% confidence level. The ratio found for the spectrum 
G88.1, which has the highest flux and the softest spectrum, 
is the largest, in contrast with the expectation that it should 
be the smallest since a smaller contribution of the reflection 
would result in a softer spectrum. This is most likely due to 
the coupling of the photon index, absorption, and reflection 
ratio in the spectral fitting. 

The parameterizing spectral fit results and hardness 
ratios show that the 1989 spectrum is significantly flatter 
than that from 1988. Assuming that the data can be 
described adequately by the model in equation (1), we can 


TABLE 2 

Reflection Model Fits to Spectra 


Equivalent 




Photon 



Reflection 

Line 

Width 

xVdof 

Spectrum 

(10^> cm-=) 

Index 

N(powcr law)* 

N(reflection)* 

Ratio 

Flux" 

(eV) 

G88.tot 

1.3 < 3.1 

2.0 ±0.1 

1 4 + 0.3 

2.6!?;i 

1.9!?:? 

5.0!};! 

120 ±40 

14/28 

G88.1 

3.0! i;S 

2,1 ± 0.2 

2.0!g;| 

5.7!JV 

2.9!?:*' 

4.6!!:| 

90 ±50 

32/28 

G88.2 

0< 1-8 

1.9 ± 0.1 

I-U_o.i 

1.3!!:! 

i.2!!:S 

5.0 ± 2.4 

140 ± 70 

13/28 

G88.3 

0<4.4 

1 0 + 0.3 

i.2!S:l 

2.1!®:? 

1.8!?:? 

5.o!!:t 

130!?? 

16/28 

G89 

2.2 < 13.0 

1 ./- 0.2 

0.5 !J;? 

0<60 

0<6'= 

4.0 ± 3.0 

180 ± 140 

24/28 

Joint fit 


1.96!g;g| 


2.4!?;? 

... 

4.6 ± 1.5 


39/57 

G88.tot 

1.1 < 2.5 


... 

i.8!J:| 


no ±40 

14 

G89 

5.6!!:? 


0.64!®;°! 

... 

3.8!?:!' 

... 

210 ± 70 

25 

Joint fit 

1.7!l;f 

2.oo!g:l| 


... 

. . . 



66/85 

G88.1 

1 - 8 - 0.2 

2.7!?;! 

i-6!J:J 

4.6!!:? 

100!fg 

34 

G88.2 

» * « 

1 • • 

t*-o.i 

3.4!?:l 

2.9!*:! 

5.0!!:? 

130!?g 

16 

G88.3 

• • t 

• » 4 

1.3 ± 0.2 

3.5!?;? 

2.7!?:? 

5.1!!:? 

120 ±60 

17 


NoTES.--In the second and third panels, parameters from simultaneous fits of spectra arc pven. In the top row in each case the values of 
parameters constrained to be equal arc listed, and in the succeeding rows the values of parameters left free and corresponding to each 
spectrum arc listed. 

• Model normalization ( X 10”* photons keV' * cm”* s'* at 1 kcV). 

** Line flux ( x 10” ® photons cm” * s” *) in the line, 

* Parameter upper limits constrained by the table model limits. 
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test the hypothesis that variability in a single model param- 
eter dominates the spectral variability by fitting the 19^ 
spectrum and the average 1988 spectrum simultaneously, 
tying most parameters together, and allowing various 
parameters to vary independently in turn. Doing so shows 
that variation in the absorption rather than the reflection 
ratio or the index produces the largest reduction in x (from 
48/59 degrees of freedom [dof] to 38/55 dof). These results 
are given in the second part of Table 2, and the x^ contours 
computed for the absorption and index are shown in Figure 
4. These show that long-term spectral variability assuming 
fixed reflection ratio is present at the 90% confidence level 
and can be modeled by a constant index but a change in 
absorption. 

The hardness ratio shown in Figure la shows that the 
significant spectral variability occurring during the 1988 
observation can be characterized by a change in ratio 
between the medium-energy band dominated by the power- 
law continuum and the hard energy band in which the 
reflection component contributes substantially. The charac- 
ter of the spectral variability suggests that either a change m 
index or reflection ratio has occurred. Fitting all three 
spectra simultaneously and requiring the spectral shape to 
stay the same with only a change in total normalization to 
account for the flux change gives x^ = 87.8/87 dof. Again 
assuming that the spectral variability is dominated by 
variability of only one parameter, fitting with equation (2) 
shows that allowing either the photon index or reflection 
ratio to vary gives an equally good fit (x^ = 66/86 dof). The 
X^ contours of spectra from the first and second days are 
shown in Figure 5, illustrating that the spectral variability 
can be explained by either a variable index or variable 
reflection ratio. Variability in the absorption explains the 
spectral variability less well (x^ = 78/86 dof). Assuming it is 
the reflection ratio rather than the photon index that 
changed, fitting with the model described in equation (1) but 
allowing only the power-law normalizations to vary inde- 
pendently shows that these data are consistent with the 
reflection flux having remained constant during the 3 days 
of observation (x^ = 69/90 dof). 

The ambiguity of these results is inherent in Ginga data 
from moderate-flux Seyfert 1 galaxies, because the decrease 



Pjq 4 — contours (90% and 1 a) for the average 1988 Ginga spec- 
trum and the 1989 Ginga spectrum for a joint fit of the reflection model, 
allowing the index and absorption to vary. These show that in the context 
of the reflection model, the spectral variability can be explained by varia- 
tion in the neutral absorption column. 



Fig 5.— contours (90% and 1 a) for Ginga spectra G88.1 and G88.2 
from a joint fit of the reflection model, allowing the index and reflection 
ratio to vary. These show that in the context of the reflection model, the 
spectral variability can be attributed to variability in the ratio or the index 
or both. 

in detector effective area toward low energies means that 
the photon index in the reflection model is not well deter- 
mined. Recent ASCA results may provide a partial prece- 
dent for the assumption that the sp»ectral variability in this 
case was produced by a change in the reflection ratio rather 
than in the power-law index. During a 3 day ASCA obser- 
vation of NGC 7469 in 1993, both 30% total flux variability 
and spectral variability was observed (Guainazzi et al. 
1994), but spectral variability was confined to lower ener- 
gies inaccessible to Ginga, while the power-law index deter- 
mined above 2.5 keV remained constant. Note that large 
changes in the photon index are detectable in ASCA 
spectra, as shown in other observations of Seyfert 1 galaxies 
(e.g., Mrk 766; Leighly et al. 1996). 

3.1.2. Partial Covering Fits to Ginga Spectra 

It was shown by Piro et al. (1990) that the 1988 Ginga 
spectrum could be modeled by partial covering by a thick 
absorber. This model fit acceptably all the Ginga spectra 
separately, although for the 1989 data, the covering fraction 
could not be constrained. This model could potentially 
explain the observed variability because if the intrinsic 
emission were constant, an increase in the fraction covered 
would result in a decrease in observed emission as well as an 
overall hardening of the spectrum. Fitting the three 1988 
Ginga spectra simultaneously and allowing only the cover- 
ing fraction to vary independently gave a statistically 
acceptable fit (x^ = 78/87 dof) with the covering fraction 
varying from 0.098 during the first day to 0.37--0.32 during 
the second and third days. However, determining the con- 
tribution to the total x^ from each spectrum separately 
shows that this fit was unacceptable for the spectrum from 
the first day (x, = 1.56/27 dof). Allowing the partial cover- 
ing column density to vary also did not result in a signifi- 
cant improvement in fit. Thus, if partial covering by neutral 
material is the cause of the flux and spectral variability, the 
intrinsic luminosity is required to vary in concert with the 
covering fraction, an unlikely scenario. 

3.2. The ROS AT Data 

The hardness ratios shown in Figure 2 and the PHA ratio 
shown in Figure 3 demonstrate that spectral variability 
occurred between the two ROS AT observations, but not 
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during each observation. Accumulating the data from each 
observation and fitting with a single absorbed power law 
from 0. 1-2.0 keV reveals evidence for spectral complexity 
(R91: Xy = 1.41/27 dof; R92: xl = 1-61/27 dof; see also 
Turner et al. 1993a, Brandt et al. 1993). Addition of an 
absorption edge, shown to be generally appropriate for 
ROSAT spectra from Seyfert 1 galaxies (e.g.. Turner et al. 
1993a), improves the fit substantially (Table 1). The column 
densities were consistent with each other and were slightly 
larger than the Galactic value of 4.8 x 10^® cm”^ (Elvis, 
Lockman, & Wilkes 1989). The photon indices were signifi- 
cantly different, much larger than those found in the Ginga 
spectra, and they were correlated with flux as suggested by 
the softness ratios in Figures 2a and 2b. The absorption 
edge energies, consistent with absorption by highly ionized 
oxygen, and optical depths were not well determined and 
are consistent with each other. 

It is well known that some time-dependent residual 
uncertainty in the PSPC calibration remains. Although the 
gain change is a nonlinear function of energy, the change in 
gain seems to be a linear trend with time {ROSAT bulletin 
64). These observations were made just 6 months apart, 
suggesting that the gain change between them should not be 
large. Further, the systematic error should be no larger than 
the difference between the two currently released ROSAT 
PSPC-B (Position Sensitive Proportional Counter) detector 
response matrices (referred to as “gainl” and “gain2”). 
Thus, the effect of a possible gain change can be tested by 
fitting each spectrum using the two different matrices simul- 
taneously. Doing so using the above model and computing 
the contours shows that model parameters are consistent 
for the two matrices but still differ significantly between the 
two spectra. Thus, the spectral variability observed between 
the two ROSAT spectra cannot be attributed to the residual 
uncertainty in the PSPC calibration. 

3.2.1. Soft Excess Modeling 

When a blackbody is used to model the spectral complex- 
ity (Table 3), the resulting power-law indices of 1. 8-2.0 are 
consistent with the mean for Seyfert 1 galaxies found after 
reflection is accounted for (Nandra & Pounds 1994); in 
contrast, this model could not explain the ROSAT spec- 
trum from NGC 3783, as the resulting extremely flat index 
was considered physically implausible (Turner et al. 1993b). 
The absorption becomes pegged at the Galactic value in the 
spectral fitting of the 1992 data, possibly an indication that 
the soft excess component is broader than a single black- 
body (e.g., NGC 5548; Done et al. 1995). Fitting with an 
additional blackbody results in a reduction in x^ by 7.3 with 
a reasonable absorption column {N^ = 5.4 x 10^^ cm“^) 
but a flatter index which could no longer be constrained. 
The absorption, photon index, and blackbody temperature 
are consistent between the two ROSAT spectra. Fitting 
jointly and leaving only the continuum normalizations free 


gave a good fit (x^ = 50/50 dof) suggesting that most of 
observed spectral variability can be explained by a cha 
in relative normalization between the power law and 
blackbody component. Allowing any of the three o- 
parameters to vary independently resulted in a bettei 
~ 35/49 dof for all three), but because of the p 
energy resolution and limited bandpass of the PSPC 
cannot be determined which parameter further contrib 
to the spectral variability. 

3.2.2. W arm Absorber Model Fits 

The ROSAT spectra can be successfully fit using a me 
which includes an edge presumably resulting from absi 
tion by partially ionized material in the line of sight, 
ionization state of the material can be described by 
ionization parameter U = Q/Anr^cN^y where Q is the flu 
ionizing photons, r is the distance from the source of i 
izing radiation, and is the density of the ionized 
(Yaqoob, Warwick, & Pounds 1989). The photonized gn 
the line of sight should respond to continuum flux chai 
resulting in spectral variability (e.g., Halpern 1984). Ge 
ally speaking, due to reduced soft X-ray opacity ROS 
spectra should be observed to soften as the flux increa 
and the measured edge energy may increase correlated ^ 
the dominant ionization state of oxygen. In Ginga dati 
which the lower end of the bandpass is significantly at 
the oxygen edge energies, the primary result is an appa: 
correlation between absorption and flux (e.g.. Nan 
Pounds, & Stewart 1990). 

Spectral fitting was done using a warm absorber t. 
model (Yaqoob et al. 1989). Fitting the two ROSAT spe 
separately showed that the photon indices are consis 
while the log (IT) and log (N,^arJ values vary margin 
such that the ionization parameter is correlated with 
flux. The results are given in Table 4. Assuming that 
index is constant allows tighter constraints on remaii 
parameters, showing that the ionization parameter 
ionized column increase while the neutral absorpt 
decreases as the flux increases. Note that there is still str 
coupling of the three remaining parameters, and given 
low resolution of the PSPC spectra, these results are hij 
model dependent. 

3.3. The Ionized Disk Model 

The hard X-ray tail observed in Ginga spectra is thoi 
to result from scattering of incident primary X-rays b} 
accretion disk. If the reflection region is close to the nucl 
the disk surface material may become ionized, resultin. 
reflection at soft X-ray energies due to reduced opacit 
soft X-rays. Emission lines should also be produced thro 
recombination. This model without accounting for 
emission was used to explain the X-ray spectrum of Cyg 
by Done et al. (1993). Ionized disk models accounting 
line emission have been developed by Ross & Fabian (P 
and Zycki et al. (1994). 


TABLE 3 


Blackbody Model Fit Results 


Spectrum 

Nh 

(10^‘ cm-*) 

Photon 

Index 

N(power law)* 

(keV) 

iV(bIackbody)** 

zVdof 

R91 

R92 

0.5 < 0.6 
Galactic < 0.53 

2-Olt-* 

\ 17 + 0,13 
'- 0.06 

0.70!«;g| 

010!S;gJ 
0.11 ±0.01 

1.1 ± 0.2 

mm 


■ Power-law normalization ( x 10"^ photons keV" * cm”* s' * at 1 kcV). 
^ Blackbody component normalization ( x 10”’* ergs cm”* s' *). 
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TABLE 4 

Warm Absorber Model Fit Results 




Photon 




Spectrum 

(10^' cm-^) 

Index 

logc; 


zVdof 

Index free: 






R91 

• 0.531J:" 

2.2!“:} 

0.4 < 0.6 

22.7 < 23.2 

23/24 

R92 

■ 0.59ig;2’ 

2.2 ± 0.1 

0.2!“:i 

22.3!S:J 

9/24 

Index constrained: 






R91 

• 0.51!“:" 

*•1-0.1 

0.4 ± 0.1 

22.9!g;’ 

23/25 

R92 

0.58 ± 0.03 


o.2!»:i 

22.4!“:g 

10/25 


TABLE 5 


Ionized Disk Fit Results 




Photon 




Spectrum 

(10=‘ cm-^) 

Index 

f 

Ratio 

zVdof 

GSS.tot 


7 ac + 0.05 

*.o;)_o.o8 

23 < 137 

1.7 ± 0.7 

13/24 

R91 

066!J:>| 

X3 ± 0.2 

lOOlJg 

2-1!}:? 

23/24 

R92 

0.68!ggJ 

2.1 ±0.1 

100!}“ 

2-0!J:| 

14/24 


We use the model developed by Zycki et al. (1994) for 
spectral fitting parameterized by the ionization parameter 
i = 4nF/n^ where F is the integrated ionizing flux and n is 
the number density. Results are given in Table 5. The 1988 
Ginga spectrum could not constrain the ionization param- 
eter, so it was consistent with reflection from neutral 
material. The ionization parameter could be constrained by 
the ROSAT spectra, primarily because of the spectral com- 


TABLE 6 

Joint ROSAT and Ginga Spectral Fit Results 
A. Warm Absorber 


Spectrum 

(10"‘ cm”") 

Photon 

Index 

Reflection 

Ratio 

Line 

Flux* 

EW 

(eV) 

log U 

'Og ^..nn 

zVdof 

Warm absorber: 

Joint 

GSS.tot 

R91 

0.53 ± 0.03 

*5 17 + 0.06 
' - 0.05 


i9!i:“ 

60 ±40 

0.4!“:J 

22.6!“:? 

41/49 

17 

24 

Joint 

GSS.tot 

R92 

0.58 ± 0.02 

2,13 ±0.04 

1 7 +0.9 

•5-5 -0.8 

3.9!};? 

90!}“ 

0.2 ± 0.1 

22.3 ± 0.2 

24/49 

14 

10 


B. Ionized Disk 


Spectrum 

(10^‘ cm-^) 

Photon 

Index 

Reflection 

Ratio 


zVdof 

Joint 

GSS.tot 

R91 

0.64!“:“! 

1 07 + 0.06 
^•^'-0.04 

l-2!gj 
4-6! 11 

120!}“ 

48/50 

15 

33 

Joint... 

GSS.tot 

R92 

0.67 ± 0.04 

2.00 ± 0.04 

2 7q??.05 
*•'^-0.69 

no ±40 

33/50 

15 

IS 


C. Blackbody 


Spectrum 


Nu 


(10^1 cm-") 


Joint Galactic <0.5 

GSS.tot 

R91 


Joint 

GSS.tot . 
R92 


0.5l!g;g| 


Photon 

Index 


1.91 ±0.04 


Reflection 

Ratio 


1 - 412:1 


Line 

Flux* 


EW 

(eV) 


4.9 ± 1.7 120 ± 40 


4.9 ± 1.7 120 ± 40 


kT 


A 1A + 0.004 


0.11 ±0.01 


Normalization** zV^of 

38/53 

15 

2-51S1 23 


0.9 ± 0.2 


26/53 

15 

II 


second and third panels, parameters from simultaneous fits to spectra are given. In the top row in each case the values of 
parameters constrained to be equal are listed, and in the succeeding rows the values of parameters left free and corresp7ndinV?o w 
spi^trun. are hst^ Note that the iron line is already included in the ionized disk model so no separate line flux is given ^ 
•Lmeflux(xl0 ’photonscm-Js-')intheline. 

the sJirS^lihs where L„ is the source luminosity in units of 10” ergs s' >. and is the distance to 
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Fig. 6a 



P”'"" law-blackbody-reflection-Unc model to the nonsimultaneous ROSAT and Gi^ga spectra. 


plexity. The resulting photon index for the 1992 spectrum 
was consistent with the Ginga result but was significantly 
larger for the 1991 spectrum. The ratio of the reflected to the 
incident flux indicates that more than half the incident flux 
must be reprocessed by the ionized material, 

3.4, Joint Model Fits 

The complexity of the broadband spectrum makes devel- 
opment of a coherent model for the emission and variability 
of NGC 7469 difficult when spectra from each mission 
separately cover only a limited energy band. An attempt is 
made to develop a coherent model by fitting together the 
nonsimultaneous ROSAT and Ginga spectra to the models 
investigated above. This procedure is justified if we assume 
that the photon index does not change and that we have 
modeled correctly the spectral variability in the Ginga data 
with a variable reflection ratio. These results are given in 
Table 6; in the first row of each group, the values from 
parameters tied together are given along with the total 
while subsequent rows display the parameters specific to 
each spectrum and the contribution from that spectrum to 
the total X • Since all the models gave good fits, they must 
be distinguished by whether the fit parameters are physi- 
cally reasonable. 

Fitting the spectra with a power law plus reflection and 
warn absorber model gives the results in the first section of 
Table 6. Reasonable values of log {U) and log con- 

strained pnmarily by the ROSAT spectra are found; 
however, the ROSAT spectra force a steep photon index to 
nt the warm absorber, which is compensated for in the 
Ginga spectrum by an unphysically large reflection normal- 
ization ratio. Further, the intrinsic index is required to be 
^ .2, which is rather more steep than that generally found 
rom Seyfert 1 galaxies (e.g., Nandra & Pounds 1994), Thus, 
e warm absorber alone does not describe the complex soft 
-ray spectrum of NGC 7469 well. Note that no evidence 
or a warm absorber was found in the higher resolution 
spectrum from the 1993 ASCA observation of NGC 7469 
(Guainazzi et al. 1994). 

ionized disk model obtained for 
e ROSAT and Ginga spectra when fitted separately as 


shown in Table 4 are mostly consistent with one another 
except for the ratio of reflected emission. Fitting the spectra 
together and allowing the ratio to vary independently gives 
consistent results; however, the fits to the ROSAT spectra 
are somewhat worse than for other combined model fits. 
This is because the ionized reflection spectrum is dominated 
by soft X-ray lines rather than the steep soft X-ray contin- 
uum required to fit the 7 spectra. 

ui soft component using a single-temperature 

blackbody produced a good fit to the ROSAT and Ginga 
sj^ctra together. In both cases, the photon index was about 
1.9, reasonable for Seyfert 1 galaxies (e.g., Nandra & 
Pounds 1994), and the column density was near the Galac- 
tic value. The reflection ratios were slightly higher than 1 
but this may be expected since the average 1988 Ginga spec- 
trum was used. The line equivalent width of 120 cV is 
consistent with production in an accretion disk or other 
cold extended material (e.g., George & Fabian 1991 ; Sivron 
& Tsuruta 1993; Krolik, Madau, & Zycki 1994). The black- 
body temperature of '^110 eV was consistent between the 
two ROS.4 7 spectra, while the flux varied by nearly a factor 
of 3 in this parameterization. Uncorrelated variability 
between the soft excess and power-law continua normal- 
izations was also found in the 1993 ASCA observation 
^uainazzi et al. 1994). The residuals and model shown in 
Figure 6 arc given for the 1992 ROSAT and the 1988 Ginga 
spectra. The soft excess could also be modeled equally well 
by a multiple blackbody accretion disk spectrum, requiring 
a slightly higher column density and temperature, but the 
uncertainties in the normalization for this parameterization 
change in flux between the two 
7 observations, showing that soft excess variability is 
model dependent. 

4. DISCUSSION 

4.1. Constraints on the Origin of the Hard X-Ray Reflection 

The significant spectral variability found in the 1988 
Gmga observation could be attributed to either a variable 
photon index or a variable ratio of reflected to primary 
emission. If the photon index is assumed not to have varied 
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over the 3 day period, then the reflection ratio must have 
changed in order to produce the spectral variability. The 
reflection flux is consistent with a constant value; the 
assurnption that it was constant leads to constraints on the 
location of the reflection region. 

One attractive explanation for the observed spectral 
variability is that much of the reflecting material is located 
at least 3 light days from the source, and so that contribu- 
tion can be determined. Two different scenarios must be 
considered, since the reflection ratio measured on the first 
day was larger than 1. That value of the reflection ratio 
might correspond to the physical reflection ratio, implying 
that the reflecting material subtends a correspondingly 
larger solid angle than 2 tc sr. Given the change in intrinsic 
emission, a constant reflection flux would imply a reflection 
ratio for the second and third day of about 0.6, implying 
that at least 40% of the reflected emission must come from 
farther than 3 light days from the center. On the other hand. 
It could be that the physical reflection ratio is 1, implying 
that the intrinsic X-ray emission was larger before this 
observation. By the same reasoning, more than 60% should 
e coming from a region more than three light days from 
the central region. 

Matt, Perola, & Piro (1991) show that for a source 
located 20r, above the disk, for inner and outer disk radii of 
u lOOOr,, respectively, and moderate inclination, half 
the reflected photons must come from the inner ~50r,. For 
a lag of 3 light days, these limits imply a black hole mass of 
- 5 X 10 M©. If the source is confined within 6r„ the light 
crossing time is about 30,000 s. The light curve shown in 
Figure Id and also the light curve shown in Figure 1 of Barr 
(1986) indicate that the doubling timescale is about 7 hr, 
just within this limit. Thus, the variability timescales indi- 
cate that it is j'ust possible that the observed constant reflec- 
ted emission could have come from the relativistic disk. An 
alternative site for the reflection is the molecular torus, 
hypothesized to lie outside the broad-line region. Recent 
modeling by Ghisellini, Haardt, & Matt (1994) and Krolik 
demonstrates that an iron fluorescence line and 
reflection hump can be produced, very similar to those 
found from a source located above a flat accretion disk the 
geometry considered by Matt et al. (1991) and others. 

The profile of the iron Ka line, which should also be 
produced in the material responsible for the reflection, can 
be used to distinguish between these two origins. Both the 
relativistic accretion disk and material further away such as 
a molecular torus are expected to produce a line of ~ 100 
eV equivalent width when illuminated by a hard X-ray con- 
tinuum (e.g., George & Fabian 1991; Krolik et al. 1994). If 
the iron line is produced in the relativistic accretion disk the 
profile should be broadened by gravitational and Doppler 
effects (e.g.. Matt et al. 1991), while if the line is produced in 
the molecular torus, the line should be narrow. The ASCA 
spectrum of NGC 7469 had a narrow iron Ka line with 
(Guainazzi et al. 1994). This contrasts with the 
II of Seyfert 1 nuclei, which contain 

clearly broad lines (e.g., Mushotzky et al. 1995). In the long 
observation of MCG —6-30-15, the line profile clearly sug- 
o emi^ion from a relativistic accretion disk (Tanaka et 
al. 1995). Thus, the observed narrow line in NGC 7469 
provides further evidence that the reflection comes from a 
more extended re^on. An equivalent width change corre- 
sponding to the reflection ratio change is also expected, but 
It cannot be constrained with these data. 


The lack of change in the reflection component flux was 
suggested by observed spectral variability in Gittga data 
from two other Seyfert 1 nuclei: Mrk 841 (George et al. 
1993) and NGC 5548 (Nandra et al. 1991). In both these 
cases, though, the interval between the observations was far 
too long to place any useful constraint on the geometry. 


4.2. TheT emperature of the Thermal Component 

The soft excess component in the ROSAT data was 
described most successfully by thermal models. The tem- 
peratures found ranged between 100-110 eV for a single 
blackbody and ~ 140 eV for the inner edge temperature of 
the multiple blackbody accretion disk spectrum. The black- 
body temperature for the soft excess component in the 
ASCA spectrum was found to be 80(66, 99) eV (Guainazzi et 
al. 1994). These temperatures may be slightly high to be the 
ernission from a bare accretion disk (e.g., Ross, Fabian, & 
Mmeshige 1992), but they are consistent with thUe 
observed in ROSAT spectra from other Seyfert 1 nuclei 
with similar luminosity (e.g., Walter & Fink 1993). 

When the ROSAT and Ginga data were fitted simulta- 
neously with a power law, reflection, and blackbody model, 
it was found that the spectral variability in the ROSAT data' 
coujd be explained by a change in the relative normal- 
ization between the blackbody component and the power 
law, so that the softer brighter ROSAT spectrum could be 
explained by the greater relative contribution of the soft 
excess component. In the ASCA observation, the spectrum 
was found to soften with an increase in flux, implying that 
the variability of the two components was not correlated 
(Guainazzi et al. 1994). Two-component spectral behavior 
was also found in the EXOSAT data, where the 1-7 keV 
X-rays varied significantly while the 0.05-2 keV energy 
remained constant (Barr 1986), which has been sub- 
sequently explained by a constant thermal component 
(Walter & Courvoisier 1992). These results may imply that 
thermal reprocessing of primary X-rays in accretion disk 
material, posited to explain correlated variability found in 
the UV-optical bands from some Seyfert 1 galaxies (e.g., 
Clavel et al. 1992), is relatively less important on long time- 
scales m this object and could be further evidence that hard 
X-ray reprocessing does not occur primarily in the inner 
disk. 


4.3. The Absorbing Material 

We found that the best explanation for the flatter spec- 
trum of the 1989 Ginga observation was a change in the 
neutral absorption. Although the statistics were poor, this 
result is interesting, since similar inverse correlations of 
absorption with flux have been seen between Ginga obser- 
vations of other Seyfert 1 galaxies, including Mrk 335 
(Turner et al. 1993c). 

Barr (1986) noted that because of the large reddening 
seen in the stellar spectrum and narrow line emission, the 
soft excess in the EXOSAT spectra may be much larger 
than measured. The reddening of E{B-V) = 0.3 would 
correspond to a neutral absorption column of Nj, = 1.8 
X 10 cm assuming a standard dust-to-gas ratio. The 
low bandpass of ROSAT allows precise measurement of the 
absorption column, especially considering the good sta- 
tistics of the spectrum from 1992. Although the model for 
the soft component is uncertain, none of the models con- 
sidered here required a neutral absorption column of more 
than 7x10 cm , less than one-third of the column esti- 
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mated from reddening. Thus, there is no evidence that gas 
corresponding to the optical reddening covers the nucleus. 

ternatively, a gas-to-dust ratio approximately 3 times 
Galactic may be implied. 

A warm absorber can explain the ROSA T spectra alone, 
but It IS less successful in explaining the combined ROSAT 
Further, in the higher energy resolution 
f j evidence for a warm absorber was 

oun ( uainazzi el al. 1994). This result contrasts with 
those from other ASCA observations of Seyfert 1 galaxies 
which find soft X-ray edges explained by absorption 
ionized material to be common (e.g., 
MCO -6-30-15; Fabian et al. 1994; Otani 1995). However, 
repeated observations using ASCA show that the column 
density of the warm absorber is sometimes variable (Fabian 

^ observation of a warm 

a sor er during a single observation does not preclude its 
existence. It is possible that the larger neutral absorption 
column discovered in the 1991 Ginga observation could be 
* ^ warm absorber (e.g., Nandra et al. 1990), 
with the statistics of the spectrum being too poor to detect 
the ionized iron absorption edge. Further, the 40% flux 
variability observed in 12 hr during the 1992 ROSAT 
° with no accompanying spectral variability is 

probably comistent with the warm absorber model, since it 
wou d be difficult to measure fine changes in the model 
parameters considering the energy resolution and statistics, 
and since the ionization state may not change on timescales 
^ than about 1 day. Other reasons why warm 
absorbers may not be observed to vary are explored in 

Turner etal. (1993b). 

5. CONCLUSIONS 

We have presented analysis of the ROSAT and Ginga 
observations of NGC 7469 in order to address the X-ray 
emporal and spectral variability. The most notable result is 
me detection of significant spectral variability in the hard 
X-rays occurnng during the 1988 Ginga observation. Since 
e spectrum can be most successfully explained by a power 
aw plus reflection model, to produce the observed spectral 
vanability a change in either the power-law index or the 
retlection ratio must hav6 occurred. Assuming that the 
power-law index did not change (an assumption based par- 
la y on the observed lack of change during an ASCA 
observation; Guainazzi et al. 1994), the reflection ratio must 
nave changed. The reflection flux could not be well deter- 
mined but is consistent with a constant value. This fact, plus 

Ka line observed in the 
ASCA .spectrum (Guainazzi et al. 1994) makes an extended 


origin for the reflection, for example, in the molecular i 
an attractive hypothesis. 

Spectral variability was found between two RC 
observations 6 months apart, which was clearly much 
from uncalibrated gain chanj 
me UbFC. The spectral complexity could be modeic 
various soft components, including a blackbody a 
warm absorber. Fitting nonsimultaneous ROSAT 
Ginga spectra together, again under the assumption th; 
photon index had not changed, found the black 
thermal component favored as the explanation for the 
excess. This explanation is supported further by the o 
vation of a soft excess component and no ionized ox 

® spectrum (Guainazzi\ 

jh observed spectral variability could thei 
exp amed by a change in the relative normalization bet' 
the power-law and thermal components. 

• ^ r^ommon assumption in this analysis is that the ph 
index ‘!*‘f"ot change. This assumption is necessary bee 
ol the limited bandpass of both the PSPC and the LAC 
ambiguity of determining the photon index in Ginga sm 
Will also be a problem in the analysis of X-Ray Tir 
Explorer (XTE) spectra from Seyfert 1 objects exhib 
‘I'is can be overcome successful!' 
coordinated observations with ASCA because of' 
mcmase m effective area toward low energies providec 
YTC monitoring experiments using Al 

tkr« ^ Sood way to determine dirt 

through variability whether the relativistic disk produ 
the broad line m other Seyfert 1 galaxies is also the origi 
the reflection spectrum. ^ 
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